Clock shifts in a Fermi gas interacting with a minority component: a soluble model 
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We consider the absorption spectrum of a Fermi gas mixed with a minority species when majority 
fermions are transferred to another internal state by an external probe. In the limit when the 
minority species is much more massive than the majority one, we show that the minority species 
may be treated as static impurities and the problem can be solved in closed form. The analytical 
results bring out the importance of vertex corrections, which change qualitatively the nature of 
the absorption spectrum. It is demonstrated that large line shifts are not associated with resonant 
interactions in general. We also show that the commonly used ladder approximation fails when the 
majority component is degenerate for large mass ratios between the minority and majority species 
and that bubble diagrams, which correspond to the creation of many particle-hole pairs, must be 
taken into account. We carry out detailed numerical calculations, which confirm the analytical 
insights and we point out the connection to shadowing phenomena in nuclear physics. 
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I. INTRODUCTION 

The spectroscopy of spin excitations in atomic systems 
is important for basic science as well as being technolog- 
ically relevant to atomic clocks. The subject has a long 
history, going back to studies of spin-exchange optical 
pumping [l[ and of line shifts in hydrogen masers Q . In 
recent years it has acquired renewed interest following ex- 
periments on ultracold atomic gases that have played an 
important role in probing effects of interatomic interac- 
tions in these systems [3|— 16| . In a typical experiment, one 
induces transitions of atoms from one hyperfme state of 
the ground state manifold, denoted by 1, to a second hy- 
perfine state, 2, in the presence of atoms in a third state, 
3. Particular interest has focused on situations where the 
interatomic interactions are strong, for example for 6 Li 
for which scattering lengths have magnitudes ~ 10 3 ao for 
a large range of magnetic fields. 

The quantity measured in experiment is basically a 
two particle correlation function that is difficult to cal- 
culate when interactions are strong. Many effects have 
to be considered, including particle self energies, vertex 
corrections, pairing, and the inhomogeneity of the atomic 
cloud LZHia. In this paper we consider a simple model 
where the mass of the bystander atom, 3, is much larger 
that of states 1 and 2. This allows us to include self en- 
ergy and vertex corrections to all orders in a conserving 
approximation which becomes exact when the system is 
highly polarized, in the sense that the density of the by- 
stander atoms is much smaller than the density of the 1 
atoms. Throughout most of the paper we shall neglect 
the interaction between 1- and 2-atoms which does not 
give rise to shifts in the absence of interactions with 3- 
atoms. For densities of 3-atoms low enough that they are 
nondegenerate, the statistics of these atoms plays no role, 
so our calculations apply to both bosons and fermions. 



For the case of nondegenerate majority atoms, our for- 
malism enables us to derive in a straightforward way ana- 
lytical results obtained previously 0, Q] . The calculations 
bring out the important role of the processes analogous to 
those considered by Aslamazov and Larkin fl7j in studies 
of fluctuation contributions to response function close to 
the transition temperature in superconducting metals. 

We find that vertex corrections can qualitatively 
change the clock shift compared with the prediction with- 
out vertex corrections. We also show that large line shifts 
are not associated with resonant interactions. For in- 
stance, when one interaction, e.g. 1-3, is on resonance, 
the clock shift has the same magnitude but the opposite 
sign compared with its value when the 1-3 interaction is 
zero. Our analytical results are confirmed by numerical 
calculations. Another conclusion of the work is that for 
the case of massive bystander atoms, the common ap- 
proximation of including only ladder diagrams is inade- 
quate, since particle-hole correlations must be considered 
on the same footing at particle-particle and hole-hole cor- 
relations. We also discuss the relationship of the physics 
of the clock shift problem to the phenomenon of "shad- 
owing" in nuclear physics, the fact that, e.g., the total 
cross section for scattering of a pion from the deuteron 
is not equal to the sum of the cross section for scatter- 
ing from a proton and that for scattering from a neutron 

This paper is organized as follows. In Sec. II we de- 
scribe the basic formalism for calculating the transition 
rate, and in Sec. Ill we describe the calculation of the 
line shape under the assumption that the massive atoms 
may be treated as static impurities. After deriving ana- 
lytical result we present results of numerical calculations. 
Section IV is devoted to showing from diagrammatic per- 
turbation theory that for a mobile minority species with 
a large mass, the problem reduces to that of scattering 
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FIG. 1: (a) The propagator G for the fermions scattering on 
impurities, (b) The correlation function T>(u). (c) The vertex 
function. Thick solid lines indicate G, thin solid lines Go, and 
dashed lines scattering on a impurity marked by •. 



from static impurities. There we also consider the rela- 
tionship of our calculations to the X-ray edge problem 
and the phenomenon of "shadowing" in nuclear physics. 
Finally, Sec. V contains concluding remarks. 



II. TRANSITION RATE 

We consider a gas of fermions in an internal state 1 
with density n\ and mass m which interacts with a gas 
of fermions or bosons of mass TO3 and density which 
is assumed to be much smaller than ti\. The gas is sub- 
jected to a homogeneous probe that flips the fermions 
from state 1 to state 2 at a rate which within linear re- 
sponse theory is proportional to 




(1) 

where initial states are denoted by i and final ones by 
/, and their energies by Ei and Ef. The frequency of 
the applied field is 10. (We put ft and the Boltzmann 
constant equal to unity throughout.) The probability of 
occupation of the initial (final) state is denoted by Pj 
(Pf). The operator V4-( r ) creates a fermion in state a at 
position r. In terms of correlation functions, the rate is 
proportional to 

ImX>(w) oc J drdr'ImV{r,r\oj) (2) 

where T>{v, r', lo) is the Fourier transform of the quantity 
-i6(t-t')([4(v, i)Vi(r, t), ^4(r', i')V 2 (r', t')}) which may 
be regarded as the correlation function for the pscudospin 
operator that describes atoms in the states 1 and 2. 



III. STATIC IMPURITIES 

In this section, we consider the case when ms > m so 
that the 3-atoms may be treated as static impurities, as 
we shall demonstrate in Sec. IV. We shall work at nonzero 
temperature T, in which case the frequencies arc to be 
regarded initially as Matsubara frequencies, odd multi- 
ples of ttT for fermions and even multiples for bosons. 
The real time correlation function is then obtained in the 
standard way by analytically continuing from the imagi- 
nary time domain. The fermion propagators in the pres- 
ence of the impurities are 

G^zy 1 =G^,z)- 1 -£ ff (z), (3) 

where 

GKp^z)- 1 =z- P 2 /2m-e a + n a (4) 

is the bare propagator and T, a (z) the self energy. Here 
e a is the energy of a nonintcracting a- fermion (a = 1,2) 
at rest and the chemical potential. In Eq. (J3j> it is un- 
derstood that the propagator is averaged over a random 
distribution of impurities (l9| . but we shall not indicate 
this explicitly in the notation. A similar remark applies 
to the correlation function T>(ui). 

To lowest order in n 3 , the self energy has the form 

Y> a {z) = n 3 %(z). (5) 

Here T a (z) is the T- matrix for scattering of a a- fermion 
on an impurity, which is given in a matrix notation by 

%(z) = V a + V a G a (z)TAz), (6) 

where the momentum sums are implicit. We have as- 
sumed that the range of the interaction V a between the 
impurities and the fermions is much shorter than the 
lesser of the typical interparticlc distance and the ther- 
mal de Broglie wavelength, (2-K/mT) 1 ^ 2 . In this case, 
for the momenta of interest, the scattering amplitude de- 
pends only on the energy. Equations d3j)- (JSJ) are shown 
diagrammatically in Fig. HJa). 

The correlation function T>{uj) is shown diagrammati- 
cally in Fig. [IJb), and the vertex function in Fig. [TJ(c). 
The importance of vertex corrections may be illustrated 
by considering the case when the interaction between an 
impurity and a fermion is the same for the two fermion 
species. The Hamiltonian is then SU(2) symmetric with 
respect to rotations between the states 1 and 2 and the 
correlation function T>{uj) is unaffected by interactions 
Q. However, the self energy corrections to the single- 
particle propagator are nonzero and they must therefore 
be canceled by the vertex corrections. To recover the 
SU(2) invariance and to satisfy conservation laws, it is 
necessary when calculating the correlation function to 
use as the vertex in the particle-hole channel the quan- 
tity 5T,/SG (£ and G are matrices in 1-2 space) ||[2(|. 
The structure of the vertex corrections is shown in Fig. 
[IJc). We now perform such a conserving calculation of 
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T>(lo) which takes the effects of the impurities into ac- 
count exactly to lowest order in n 3 . 

The vertex corrections correspond to processes in 
which a 2-particle and a 1-hole scatter from the same 
impurity. The resulting effective interaction [the "bow- 
tie" part of the diagram in Fig. Q3c)] is given by 



V e s = n 3 T\{iio v )T 2 {i^ v + iuij), 



(7) 



since for a static impurity the energy transfer between 
particles and the impurity is zero. Because the scatter- 
ing on an impurity is independent of momentum, the 
inclusion of the vertex corrections simplifies significantly, 
and the diagrams for the correlation function may be 
summed. The result is 



(2n) 3 / d 3 p Gi (p, iuj w )G 2 (p, ioj v + iuj 1 



nzTxiiwvYFiiiwv + iw 7 )(27r) 3 J d 3 pGi(p,iuj„)G 2 (p,iuj 1J + iw 7 ) 



(8) 



r 



with ujy being a fermion Matsubara frequency and w 7 a 
boson one. The momentum integral yields 



M(zi,z 2 ) 



d 3 p 



Gi(p, zi)G 2 (p, z 2 ) 



(2ir) 3 

d 2 (z 2 )sgn(imz 2 ) - di(zi)sgn(Imzi) 
z 2 - Zi + H2 - Mi - A + Ei (21) - E 2 (z 2 ) 



(9) 



with d a {z) = m 3 ! 2 ^/ z + [i„ — e CT — E CT (z) / 'V2ir 2 and 
A = e 2 — ei the hyperfine splitting between the two 
fcrmionic states. For E = 0, d is the free particle density 
of states. We evaluate the sum over Matsubara frequen- 
cies in ((5J) by converting it to a contour integration in 
the usual way by multiplying the integrand by the Fermi 
function f(zi) = [exp(/3zi) + and choosing a contour 
that encircles the poles of the Fermi function. The inte- 
gration contour may be deformed to lie above and below 



J 



ImD(w) 



de 



(/ 2 -/i)Im 



rfi 



the cuts of the functions d\ and d 2 , which are located at 
Z\ = e and z\ — e — iw 7 , respectively, where e is real. 
After the analytic continuation — > u> + irj, with the 
physical frequency given by lj = u> + fi2 — Mij we obtain 



f°° de 

J ^-fi^lSit + vht + u + ni)- 



S(e — ir), e + Q + irj) + S(e — u — it], e + ir/) 

-S(e-Q-iri,e-iri)}, (10) 



where 



S(zi,z 2 ) 



M(z 1 ,z 2 ) 



(11) 



1 - n 3 Ti{zi)T 2 (z 2 )M(z ll z 2 )' 
The imaginary part of the correlation function becomes 



- A + n 3 [Ti -T 2 - iTrTiT 2 (d 2 - di)] cu - A + ? l3 [T 1 * - 75 - ™T?T 2 (d 2 + d*)] 



(12) 



In dUD, /1 = /(e), / a = f(e + w) and d l and 71 arc 
evaluated at the energy e+irj and d 2 and T 2 at the energy 
e + uj + i?;. The T-matrix, given by Eq. has the form 



1 + indaTai 



(13) 



where 7^ V ac is the T-matrix for a-fermions scattering at 
zero energy in a vacuum. From this it follows that Eq. 
([T2"j) reproduces the unshiftcd ideal gas result when the 
interaction is SU(2) symmetric with identical scattering 
between an impurity and the two fcrmionic species. Note 
that it is crucial to use full propagators to recover this 
symmetry. 



A. A simple limit 

Equation (fT2")l satisfies the conservation laws regardless 
of the magnitude of 713. We now study the interaction 
effects on V{uj) to the lowest order in n 3 and neglect all 
medium effects except the factor n 3 in front of the T- 
matrices in the denominator in (|T2"j) . The T-matrix is 
thus replaced by its value in a vacuum given by 



1 



2irdo 



(14) 



where 5 a is the scattering phase shift in a vacuum and 
do = m 3 / 2 ^fe/ '\f2-K 2 is the free-particle density of states. 
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For low energies, the phase shift is given in terms of the 
scattering length a a by tan S a = —ka a , where k = \j2mt. 

We make the variable change e + Hi — ei — > e, so that 
e is the kinetic energy of a 1-fermion. The 2-fermion has 
kinetic energy e + uj — A. In the limit of a low density of 
3-atoms, the line shifts arc small. One may then neglect 
differences between uj and A, and therefore the phase 
shifts of the two fermions are to be evaluated at the same 
kinetic energy. In total, keeping only the lowest order 
effects of n 3 in (fT2"j) yields 



lmP(t 



Im 



de- 



-d f(e + ei - Hi) 



- A - in 3 [ e 2i (*i-*0 - l] /27rd ' 

(15) 

We have written Eq. (|15p for the case where there are no 
2-fcrmions present initially so that the first Fermi func- 
tion in Eq. (|12|) is zero. This corresponds to a typical ex- 
perimental situation. For equal interaction between the 
an impurity and the initial 1- and final 2-fermions, i.e. 
for Si = <5 2 , ones sees immediately that the interaction 
effects vanish in the denominator of (|15[) and we recover 
the ideal gas result 



"1 



A 



IT] 



(16) 



In general, Im2?(u;) is the sum of Lorentzian lines, with 
the energy-dependent frequency shift 



Aw(e) = n 3 



it sin(2<5i - 2S 2 ) 



and with full width at half maximum equal to 



F(e) = n 3 



47T sin 2 (<5i — <5 2 ) 



(17) 



(18) 



Equations (fTT| - (|18p . which apply for arbitrary degree of 
degeneracy of the 1-atoms, have a form similar to those 
derived for a classical gas in Refs. [1 Eq. (82)] and 0, 
which studied the equation of motion for the density ma- 
trix. There is, however, a difference, since to obtain the 
result in these papers one must replace the term contain- 
ing the phase shifts in the denominator of Eq. (fT5)) by its 
thermal average. Since A and T are energy dependent, 
the line shape given by Eq. (|15p will not be Lorentzian 
in general. The magnitude of deviations from Lorentzian 
behavior will depend on the variation of [e 2 ^ 1 ~ S2 *> — l]/fc, 
over the distribution of the momentum k of 1-atoms. 

For small phase shifts, 8 a ~ ka a and Eq. (fTTj) repro- 
duces the low density expression for the shift, Aw ~ 
n327r(o2 — ai)/m, the factor of two, rather than the usual 
factor of four for the case of fermions of equal mass, be- 
ing due to the fact that we have taken the 3-atoms to be 
infinitely massive. 

Equations (fTT)) - (fl"5)) clearly illustrate the importance 
of vertex corrections. When these are neglected, the cor- 
responding results are 
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FIG. 2: The transition rate as a function of frequency and 
interaction. The inset shows the transition rate for varying 
impurity concentration. 



and 



r(e) 



n 3 - 



47T sin Si + sin S 2 



(20) 



There are a number of important conclusions that may 
be drawn from the above results. First, without vertex 
corrections the shift and damping do not display the re- 
quired SU(2) symmetry for Si = S 2 - Second, for small 
phase shifts the line width is proportional to (oi — a 2 ) 2 , 
whereas without vertex corrections the corresponding re- 
sult is of +a|. Thus, even in the limit of small phase shifts 
it is important to include vertex corrections, which give 
an interference term — 2ai0 2 . Third, the largest shifts are 
obtained for Si — <5 2 f=s 7r/4 + vn /2, where v is an integer. 
Thus, resonant scattering is not particularly favorable for 
producing large shifts. For example, take a typical ex- 
perimental situation where scattering of the fermions in 
the initial state 1 with 3-atoms is resonant while that of 
final state 2 fermions is not: the shift is then equal in 
magnitude but of the opposite sign compared with what 
it would be in the absence of 1-3 scattering, and therefore 
the magnitude of the shift is determined completely by 
the non-resonant 2-3 interaction. Finally, large widths 
and large shifts do not go hand in hand, since the largest 
widths occur when Si — <5 2 is an odd multiple of 7r/2. 



B. Numerical results 

In Figs. [UO we present numerical results for InxD(cj) 
obtained from Eq. (|12l) . The propagators used in this 
calculation are determined fully self-consistently. The 



frequency unit is Tm^/mkp with ni 
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(|18p]. and lmD{oS) plotted in units of mfci?. In Fig. [21 we 
show the transition rate for T = and with an impurity 
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FIG. 3: The transition rate as a function of frequency for 
varying initial state interaction. 




FIG. 4: Diagrams for the self energy of a 1-fermion. (a) 
Second order contributions, (b) Third-order contributions, 
(c) A general term, (d) A diagram with a general vertex 
correction. Solid lines are fermion propagators, dashed lines 
interactions with • indicating the vertices. 



density 77,3/771 = 0.1. The scattering length for the 2-3 
interaction is a 2 = —1/kp and the scattering length for 
the 1-3 interaction varies from a\ = to a\ = —100/kp, 
which is very close to resonance. For kpai = 0, the line 
shift and width are due solely to the self energy of the 2- 
atom. We see that when the 1-3 interaction is resonant, 
the main effect is to change the sign of line shift compared 
with the result for zero 1-3 interaction. This confirms the 
discussion in Sec. IIII A"l When kpai = fci?fl2 = —1, the 
unshifted ideal gas result is recovered, as it should be, 
and the small remaining width of the calculated signal 
is entirely due to a small imaginary part we have added 
explicitly to the frequency to facilitate the numerical cal- 
culations. For comparison, we also plot the result ob- 
tained for kpcii = kpa-2 = —1 when vertex corrections 
are not included. We see that, although the predicted 
line shift is small, the width is large and one does not re- 
cover the unshifted narrow line when vertex corrections 
are ignored. 

In the inset, we compare the line shape for n^/ni =0.1 
and n 3 /rii = 0.3, keeping kpai — and kpa,2 — —1 fixed. 
To ease comparison of the results for the two different im- 
purity concentrations, we have multiplied ImD(uj) by 77,3. 
The two curves largely overlap which illustrates that the 
line shift and width essentially scale with 713 in agree- 
ment with (jTTJ) - <jT5J) . Note that higher order medium 
effects coming from the self-consistent determination of 
the propagators give rise to the slight difference between 
the results for the two impurity concentrations. 

In Fig. |3j we plot the transition rate as a function of 
di keeping a 2 = 0. The line shift is large whereas the 
width is small for kpai = —1. This is in agreement with 
the conclusions reached in Sec. IIII Al from (fT7|) - (fT8)) since 
kpax = —1 corresponds to a phase shift of Si = ir/4. 
Likewise when the scattering is close to resonance with 
kpai = —100 corresponding to Si w 7r/2, the line shift is 



small whereas the width is large. Again, this agrees with 
the discussion in Sec. IIII Al 



IV. MOBILE, MASSIVE IMPURITIES 

In this section we show that in the limit 771/7713 0, 
the results we have employed in Section Hill for the self- 
energy and vertex corrections may be derived from dia- 
grammatic many-body theory for particles of finite mass. 
An important conclusion of this section is that the ladder 
approximation, which is commonly employed in treat- 
ing strongly interacting systems is inadequate to describe 
systems with 777,3 3> m when the 1-fermions are degen- 
erate. We begin by showing this for low-order contri- 
butions in perturbation theory, and then generalize the 
considerations to arbitrary order. For definiteness, we 
shall assume that the 3-atoms are fermions, but the cal- 
culations may easily be generalized to the case of bosons, 
the only difference being that the distribution function 
for 3-atoms must be taken to be the Bose distribution. 
To lowest order in the density of 3-atoms, the results are 
independent of the statistics of the 3-atoms. 



A. Second order 

In the Hartree approximation the self energy and ver- 
tex corrections are given by Eqs. ([5]) and (J7J) when the 
T-matrix is evaluated in the Born approximation, so the 
first term we shall consider in detail is the second-order 
term, Fig. L3ta). For m/777,3 — > it is given by 
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#q - gaaKj - />U) + (i - f^K+jj-, 

10 — (p — q) 2 /2m 



, (21) 



where V CT (q) is the bare interaction between a 3-atom 
and a er-fcrmion and / p = f(p 2 /2rrii + (4 — \Xi). When 
the 3-atoms are nondegenerate the 1 — f 3 factors may be 
replaced by unity, and one then finds 



.(2) 



d 3 q d 3 p' 



Vi(q)\ 



f 3 ,(l - f 1 )+ f 3 , f 1 

J P v Jp— q/ 1 •'p'+q^'p— q 



2 J P 



' (2tt) 3 (27r) 3l-lvlyl w-(p-q) 2 /2m 
(2^ |yi(<?)l W -(p-q) 2 /2rn " "'^ 



where 



r v i 2 c ) (p,o;) = 7I< 2 d ) (p,a;) + T b ^(p,u,) 



d 3 g 
(2tt) 3 



|Vi(?)| 5 



w - (p - q) 2 /(2H 



(23) 



is the T-matrix in a vacuum (calculated to second order 
in Vi). Here 
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IM?)I 5 



1 - f 1 



ladlP ' W; "" ' (27r) 3l ' iWI w- (p-q) 2 /(2m) 



(24) 



is the contribution to the T-matrix from ladder diagrams 
(particle-particle scattering and hole-hole scattering) and 



f 1 



p q 



^-(p-q) 2 /(2m) 



(25) 



is the contribution from particle-hole scattering. This 
calculation leads to two important conclusions. First, the 
presence of the degenerate 1-fermions affects scattering 
in the particle-particle and hole-hole channel and also 
scattering in the particle-hole channel, but the effects 
of occupancy of intermediate states cancel in the total, 
which to second order is given in terms of the T-matrix 
in vacuo. Second, even though the self energy may be 
written in the form 



£1 = TrTiGS 



3) 



(26) 



with the T-matrix calculated to second order in the lad- 
der approximation, it leads to an expression of the form 
(|5|) where the T-matrix contains both ladder and bubble 
contributions. Only for nondegenerate 1-fcrmions are the 
bubble diagrams unimportant compared with the ladder 
diagrams because /' C 1. By extending the above ar- 
guments to higher-order terms, one sees that the self en- 
ergy calculated from Eq. (|2l)|) with 71 calculated in the 
ladder approximation agrees with Eq. ([5]) only for non- 
degenerate 1-fermions. 



B. Arbitrary order 

In higher-order processes, a qualitatively new feature 
appears: to obtain the result ([5]), with the T-matrix given 
by Eq. one cannot use Eq. (f2l))) with the T-matrix 
calculated in the ladder approximation. To obtain the 
correct result for m/m 3 —> 0, the nth order contribution 
to the self energy has the form 



si" ) (^-r a ) = -(-y 1 )" f 

Jo 



cLtx . . .dr„_ 2 x 



Gi(n - t ii _ 2 )Gi(t„_2 - t„_ 3 ) . ..Gi(n - r a )x 
[Gain - T n -2)G 3 (T n -2 - t„_ 3 ) . . . G 3 (n - T- a )G 3 (r a - r b ) 
+ all r permutations]. (27) 

The momentum sums are suppressed for the moment to 
highlight the essential parts of the reasoning. We refer 
to this term as an nth order contribution, even though 
we use renormalized propagators, which may give rise 
to contributions of higher order in V\ when expressed in 
terms of bare propagators. The first term in the square 
brackets corresponds to Eq. ([2^1 with T calculated is the 
ladder approximation, and it does not give the result flSJ) 
for the self energy in the limit m 3 /m 3> 1. However, if 
one adds to it contributions corresponding to all possi- 
ble ways of attaching the interactions occurring at times 
Ti . . . t„_2 to the 3-bubble with vertices at r a and t&, one 
does indeed recover the result This procedure is in- 
dicated by the last line of ([2"7| and yields in third order 
the two diagrams depicted in Fig. B|b). The structure 
of a typical high order diagram generated in this way 
is illustrated in Fig. Htc). To show that (|27p produces 
the correct impurity result, we use the fact that the non- 
interacting 3-propagator for m/m 3 — > becomes 



f3g-(e 3 -M3)r 



for t < 
for r > 0, 



(28) 



where r < corresponds to the propagation of a 3-holc 
and r > to the propagation of a 3-particle. Self energy 
contributions to G 3 can be absorbed in the chemical po- 
tential /13 for low T and m/m 3 -C 1. The point is that 
for any value of r a , t&, t\ . . . t„_2 between and /3, only 
one term inside the square brackets in (|27p will have one 
hole propagator and n— 1 propagators for the 3-particles. 

This term will scale as ~ f 3 All other terms have at 

J 1 

least two hole propagators and will be suppressed in the 
limit of low concentration of the 3-particles. Note that 
the ladder diagram is not enough to include the leading 
order diagram in f 3 for any value of r tt , t&, t\ . . . r„_2 ; 
one has to include diagrams corresponding to all possi- 
ble ways of attaching n . . . t„_ 2 to the 3-bubble. In the 
limit m/m 3 — > 0, the momentum integrals in (|27[) decou- 
ple and the integral over the 3-holc line yields the density 
n 3 . We obtain 



.dr„_ix 

Gi(r b - r n _ 2 )Gi(r n _ a - r„_ 3 ) . . . Gifa - r a ). (29) 



& 



n 3 Vr / dr 2 
h 
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In frequency space this reads 

V$ l) (u)=n a Vi[V 1 G 1 (u)) n - 1 . (30) 

Summing all orders for S gives 

Ei(w) = na[l - ^iGx^)]- 1 ^ = n 3 Ti(w). (31) 

This agrees with ([5]) and we have shown that one recovers 
the correct impurity result for the self energy in the limit 
of ri3 small and m/m^ — > 0, when all crossed diagrams 
of the type illustrated in Fig. SJc) are included. 

The same argument applies to vertex corrections. Con- 
sider 1 and 2 fermions simultaneously scattering on a 
3-particle. A typical diagram needed to be included to 
recover the correct impurity result is shown in Fig. EJd) : 
For an nth order diagram, one has to include all possible 
ways of attaching the interactions occurring at t\ . . . T n 
to the n propagators in the 3-loop. In this way, the term 
where there is only one hole in the 3-loop is included for 
any value of the time arguments. This term scales as 
n,3 whereas all other diagrams are suppressed by higher 
powers of 713. When these diagrams are included to all 
orders, the effective interaction between a 1-fermion and 
a 2-fermion both scattering on a 3-atom becomes 

V eS = n 3 [l - V 2 G 2 (uj 2 )}- 1 V 2 [1 - VxG^x)]-^ 

= n 3 7i(wi)7i(wa). (32) 

This agrees with the impurity scattering result given by 

C. Higher loops and the X-ray edge problem 

So far we have considered diagrams in which there is a 
single fcrmion loop containing fermions in states 1 and 2. 
We now comment on the effect of including contributions 
with a higher number of loops. The problem under con- 
sideration in this paper has a number of points in com- 
mon with the X-ray edge problem, where the contribu- 
tions from terms containing many fermion loops change 
qualitatively the nature of the threshold behavior [2l], [22[ 
from a step function at the Fermi surface when a single 
fcrmion loop is included to a power law whose exponent 
depends on the phase shift for scattering of an electron in 
the conduction band from a deep hole. In the X-ray edge 
problem, conduction electrons scatter from a deep hole, 
which is present only for times between that at which the 
electron-deep-hole pair is created and that at which it is 
destroyed. The complications in the X-ray edge problem 
are due to the fact that the higher-order loop contribu- 
tions depend on the times at which the particle-hole pair 
is created and destroyed. In the problem under consid- 
eration in this paper, however, the heavy atoms in the 
state 3 are present for all time. The effect of the higher- 
order loops is simply to renormalize the propagator for a 
3-atom. The self energy of a 3-atom depends on energy 
but, within the approximation of a short-range potential 



made above, is independent of momentum. When higher 
loop contributions are included, the chemical potential of 
the impurities must be adjusted so that the number of 
impurities is equal to the required value. 



D. Analogy with "shadowing" in nuclear physics 

The result (|15l) has a simple interpretation, since it 
is equivalent to the statement that the self-energy of a 
particle-hole pair due to interaction with an impurity 
is proportional to e 2 *^ 3- * 1 ) — 1. Since the self energy 
is proportional to the T-matrix for scattering of a pair 
from an impurity, which is in turn proportional to S — 1 , 
where S is the corresponding S-matrix, this implies that 
S = e 2l ( S2 ~ Sl \ In physical terms, this says that the ex- 
tra phase acquired by the pair is the sum of the phase 
changes experienced by a particle in state 2 and a hole in 
state 1. The reason that vertex corrections, which cor- 
respond to interference terms, are so important in the 
present problem is that the external field creates a parti- 
cle and the hole at the same point in space. Thus if, say, 
the particle is close to an impurity, the hole will also be 
close to an impurity. If only self-energy corrections are 
included, this is equivalent to assuming that the particle 
and the hole are uncorrelated in space. 

Insight into the result for the line shift may be obtained 
^ 1 r^ing 2 usfe^ 1 ^g 1 i(4cn^^ n 2 ^ 

— sin 2(?2 (1 — 2 sin Si), 
(33) 

which implies that the energy shift is given by the real 
parts of the self energy of a 1-fermion and a 2-hole, mul- 
tiplied by factors 1 — 2 sin 2 5. To interpret this result, we 
observe that the total cross section for scattering of a 3- 
atom by a o- fermion is proportional to sin 2 5 a /k 2 , where 
k is the wavenumber of the atom. If one changes ones 
perspective and regards the process as the interaction of 
an impurity fermion with a particle-hole pair, this equa- 
tion implies that the amplitude of an impurity fermion 
at the position of the 2-hole is reduced by an amount 
~ sin 2 6i due to scattering from the 1-fermion, and like- 
wise for the amplitude of the impurity at the 2-hole. This 
is reminiscent of the experimental observation that the 
total cross section for scattering of pions from deuterons 
is less than the sum of the cross sections for scattering of 
a pion from a single neutron and a single proton, a phe- 
nomena referred to as "shadowing" . It reflects the fact 
that the neutron and proton in the deuteron are corre- 
lated, and therefore the pion field incident on, e.g., the 
proton is reduced by scattering from the neutron [l8l ]. 
The analogy between the two situations is not complete, 
however, since in the problem considered by Glauber the 
wavelength of the pion is small compared with the sepa- 
ration of the neutron and proton in the deuteron, while in 
the problem under investigation here the wavelength of 
the particle and the hole is large compared with their sep- 
aration, which is initially zero. As a consequence, where 
k 2 appears in the present problem, this is replaced by a 
factor cx< 1/r 2 >, the average of the inverse square of 
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the separation of the neutron and proton in the deuteron. 

V. CONCLUDING REMARKS 

In this paper we have solved a simple model for 
clock shifts for hypcrfine transitions between states of a 
fcrmionic atom, in the presence of a low density of much 
more massive atoms. The calculation shows the impor- 
tance of vertex corrections, which completely change the 
dependence of the shift and the width of the clock transi- 
tion on the scattering phase shifts. The calculations are 
valid for bystander atoms, either fermionic or bosonic, 
which are very much more massive than the majority 
fermions, and an important problem for the future is to 
study a finite mass ratio. 

Throughout, we have have neglected the interaction 
between 1-atoms and 2-atoms. When there are no by- 
standers, the transition has no shift and no width, and 
this result also holds in the presence of bystanders, pro- 
vided the 1-3 and 2-3 interactions are identical, since in 
that case the SU(2) invariance still holds. However, when 



the 1-3 and 2-3 interactions arc different, the line shift 
and width can be affected by the 1-2 interaction, which 
is an unexplored effect in the cold atomic gas context. 

The calculations indicate that experiments on clock 
shifts in mixtures of atoms with different masses would 
be useful. Since pairing correlations are suppressed when 
species have very different concentrations, these would 
enable one to obtain information about correlations in a 
state less complicated than a paired supcrfluid. 

An important theoretical result of our calculations is 
that it is generally not sufficient to include just ladder 
diagrams, since in the case considered, particle-hole cor- 
relations are necessary in order to recover the correct 
result for a low density of the minority component. 
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